Is there Ornstein-Zernike equation in the canonical ensemble? 
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ynP ' A general density-functional formalism using an extended variable-space is presented for classical 

^-^ I fluids in the canonical ensemble (CE). An exact equation is derived that plays the role of the 

. Ornstein-Zernike (OZ) equation in the grand canonical ensemble (GCE). When applied to the ideal 

^1^' gas we obtain the exact result for the total correlation function /ijv. For a homogeneous fluid with 

W , N particles the new equation only differs from OZ by 1/N and it allows to obtain an approximate 

expression for ftjv in terms of its GCE counterpart that agrees with the expansion of ftjv in powers 
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f-^ ; of i/N. 
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The Ornstein-Zernike (OZ) equation is an essential ingredient of any theory of the equilibrium microscopic structure 

of classical fluids since it provides an exact relation between correlation functions at the two-particle level. In the 

world of the integral equation theory p^pl , most integral equations are based on this relation closed with an additional 

approximate relation between the correlation functions. Closures like the Percus-Yevick approximation, the hypernet- 

ted chain approximation, or the mean spherical approximation, together with the OZ relation have led to well-known 

theories of classical fluids. In the world of density functional theory (DFT) [||-g, one makes explicit approximations 

j^ , for the free-energy functional and, after using a variational principle, one obtains the direct correlation function via 

Ch ' functional differentiation. Inserting the direct correlation function into the OZ equation allows to obtain the pair 

I , correlation function. Furthermore, as we shall see below, the OZ equation is an identity that arises naturally from 

X3 . the DFT formalism |-|. 

^ ' When using the OZ equation it is important to notice that it has been formulated for the grand canonical ensemble 

f_\ ■ (GCE) and thus care must be taken when considering it in other statistical mechanics ensembles. Of particular 

I— 'I interest is the canonical ensemble (CE), not only from a fundamental point of view but also because of its relevance 

, ■ in the analysis of systems with a fixed, finite number of particles N. The aim of this Letter is to study the possibility 

K^ of establishing an OZ relation in the CE Q. We shall see how the fixed- iV constraint precludes the direct application 

(^ ' of the GCE-OZ relation, this problem will be solved in a DFT context by using an extended variable-space in which, 

T-H [ together with the familiar inhomogeneous density, a Lagrange multiplier is considered as a variable. 

The OZ equation for a nonuniform system with fixed chemical potential fi and temperature T, i.e., in the GCE, in 
the presence of an external potential T4xt (r) , is usually expressed as the following relation between the total correlation 
function h and the direct correlation function c"' 



/i(ri,r2) = c(2)(ri,r2) + / p(r3)c(2^(ri,r3)/i(r3,r2)dr3. (1) 

^ ', where p(r) is the inhomogeneous density, and, in terms of functional derivatives, h and c^'^' are defined (fixed T,fi) by 

S ■ r'T^^ = -p(ri)p(r2)Mri,r2)-p(ri)<5(ri-r2) (2) 

vJ ■ oVoxt(r2) 



and 



/V;xt(ri)] _ (2), „^ <5(ri-r2) 



^. /?:i^;^p^ = c(^)(ri,r2)-:i^l^-^, (3) 



with P — l/fceT. From (|l|)-(§) one easily obtains 

5p{ri) SVcxtirs] 



SVcxtirs) Sp{r2) 



dr^ == (5(ri - r2) , (4) 



which expresses that, via the OZ equation, h and c*^^-* are (essentially) functional inverses. We note that this functional- 
inversion statement is only possible if there is a one-to-one relation between p(r) and V'oxt(r) (fixed T,fi), so that the 
OZ equation can be viewed as a consequence of this relation. This viewpoint is usual in DFT of classical fluids in the 



grand canonical ensemble (GCE) where the one-to-one relation between p(r) and Kxtlr) is a key theorem (see, e.g., 
1^) and the OZ equation follows as a corollary ||]. 

We now focus our attention on the canonical ensemble (CE) where the number of particles N is fixed. This means 
that the CE inhomogeneous density pN is normalized to N, i.e., 

PNir)dr^N. (5) 

Furthermore, from the usual definition of the two particle density p]^ {see, e.g., [Q) it is direct to obtain 
J p^^'{ri,r2)dr2 = (TV - l)pAr(ri), and thus, taking into accoimt that p^^'{ri,r2) = pAr(ri)pAr(r2)f /iAr(ri,r2) -I- Ij, 
the following constraint for h^ arises in the CE: 

PAr(ri)/iAr(ri,r2)dri = -1, (6) 

which, as pointed out by Ashcroft 0], in the context of electronic systems leads to the notion of exchange- correlation 
hole. 

We remark that Eq. (|6|) is a result of considering a system in the CE where the number of particles cannot 
fluctuate. An important consequence of this equation is that the OZ equation (]^ is no longer valid in the CE, since 
it is inconsistent with (||). Furthermore, since the OZ equation is associated to the one-to-one relation between p{r) 
and Kxt(r) it would seem that there is not a CE counterpart of this one-to-one relation, which would be of serious 
concern in DFT of the CE. The goal of this Letter is to clarify this point and to introduce a CE-OZ equation similar 
to(^. 

It is clear that, fixed N and T, the CE density pjv is determined by the external potential T4xt, but this potential 
is not (strictly) unique: any potential that differs from Vext by an additive constant leads to the same pN and 
thus the inhomogeneous density is not uniquely determined by the external potential. This does not contradict 
the first Hohenberg-Kohn theorem for DFT because it allows for this additive constant. However, this (strict) 
nonuniqueness means that the (Legendre) transform from the functional variable 14xt to the functional variable pN 
is not invertible and thus 5pAr(ri)/(514xt(r2) is a singular matrix, which in turn implies that Eq. (^ cannot hold for 
the pair 14xt, Pn- We want to stress that this is due to the fixed- A^ constraint which implies that — 6N/6Voxt{'>^2) = 
J 5yOAr(ri)/(5T4xt(i"2)rfri which shows that the matrix is singular Q. Our proposal for solving this difficulty is to 'bypass' 
the constraint by means of an extended variable space where instead of T4xt alone we consider the set {14xtj A''} and, 
instead of pN, the set {pn, ^n} where Xn will be defined later and identified with the Lagrange multiplier related to 
the fixed- A^ constraint. The role played by the different variables is clarified by rebuilding the main DFT results in 
this extended space. Our derivation is based on Legendre transform approach of Ref. H . 

The CE Helmholtz free energy F is a functional of the external potential 14xt and a function of the number of 
particles A' and the temperature T, i.e., F = F{T, N, [Kxt])- The density pjy and the 'chemical potential' Aat (different 
from pi) are functions of T and A^ and functionals of T4xt, and are obtained from the following derivatives: 

PNir,[V..t]) ^ ( ^T^] , (7) 



SV,xt{r) 




N 



AiviKxt] ={^] (8) 

where the explicit dependence on T has been been omitted. We next consider the Legendre transform of F to the 
new variables pN and Ajv, 

F{Xn,[pn]) = F{Xn,[pn]) - / PN{r)Vcxtir,XN,[pN])dr - Xn / PN{r)dr, (9) 

where we have used N — J pjv(r)dr [Eq. (0)] and we define F{Xn, [pn]) = F{J pjv(r)dr, [Vcxt{XN, [pn])])- In order to 
introduce a variational principle we 'revert' this transform by obtaining a functional that depends on both the new 
variables {pn, Xn} and the old ones {Kjxt, A^}: 



Fv.,^„NiXN, [pn]) = F{Xn, [pn]) + / PAr(r)Voxt(r)dr + AatA^ . 



(10) 



The minimization of this functional w.r.t. {pNi ^n} leads to 

fSFiXNdPN])^ 



\ SpN{r) 
and 



f dF{XN,[PN]) 

\ dXN 



+ V;xt(r) = (11) 



(12) 




Substituting (g) into ( [l^ ) and taking into account Eq. (|^) one has that the functional 

T = F + Xn pN{r)dr = F{Xn,[pn]) - / PAr(r)V'cxt(r, Aat, [pAr])dr (13) 

does not depend on Xn, that is JF = J^[pn] and thus, from ( |ll| ) we obtain 

6T[pn] 



Spir) 



+ Kxt(r) = AAr. (14) 

P=Pgc 



Comparing this equation with the GCE Euler-Lagrange equation we see that Xn plays the role of the chemical 
potential p (see, e.g., Ref. Q). If on the other hand we compare with the CE Eulcr-Lagrange equation obtained in 
Rcf . Ijl^l (see also g ) by means of the Lagrange multiplier technique we can identify Xn with the Lagrange multiplier 
associated to the constraint (^) and T is the intrinsic free energy functional in the CE. Therefore, we see that using 
the extended variable space one obtains the same Euler-Lagrange equation than in the standard DFT procedure, with 
a intrinsic free energy functional J- that only depends on the density. 

As an example we apply the above expressions to the classical ideal gas. For this system one has —(3F{T, N, [T4xt]) = 

TVlogI A~'^ /exp[— /3Vi3xt(r)]dr j — log TV! where A is the thermal wavelength. Using (|^) we obtain the familiar baro- 
metric result pAr(r, [Vi^xt]) — A^exp[— /?V;,xt(i")]// exp[— /3T4xt(r)]rfr and from (^) we would trivially obtain the corre- 
sponding result for ~PXN[Vcxt]- These results can be inverted to obtain N[pjsf] = f p]\[{r)dr and an explicit expression 
(not shown) for V^oxt in terms of pN and Ajv . Substituting these expressions in ( |l3|) we obtain the following result for 
the intrinsic free energy functional of the ideal gas in the CE: 

/3^idM = - /'pN(r)(log(A-3p^(r)) - l)dr + 0( /"pw(r)dr) . (15) 

with (/)(x) = logx! — xloga; + x. We note that this free-energy only differs from the GCE result by the term 
^(/pAr(r)dr), i.e., by (j){N). For large A^, (/'(A^) « (logA^)/2, and, since the mean square fluctuation A^N) of the 
ideal gas in the GCE is equal to A^ we find that ( |l5| ) agrees with the saddle point approximation of Ref. [|lO[ in which 
the leading correction to the CE free energy functional is given by (log[A^(A^)])/2. Finally, we recall that using (15) 
in Eqs. (P^)-(P^) we rederive the result for the equihbrium density pN, as we would have done using (|l5| ) and (14) 
with the constraint (||). 

An important property of using the extended variable space is that, as we have seen in the example, the set 
{pn,Xn} determines completely the set {V^xt,A^}, and, since the opposite is also (trivially) true, there is a one- 
to-one relation between both sets. This means that the Lagrange transform is invertible and the (Hessian) matrix 
d{pN , Xn) / d{Vcy^t, N) is not singular. Multiplying this matrix by its inverse d{Vext,N)/d{pN,XN) and equating to 
the identity matrix in the extended space, we obtain a set of equations that, in the canonical ensemble, correspond 
to the OZ equation. After some algebra we obtain the exact relations 

/iAf(ri,r2) = c^^(ri,r2)+ / pAr(r3)c^^(ri,r3)/iAr(r3,r2)dr3 -— — ^— ^ (16) 

J PN(y^2) Civ 

and 

5pjv(r2) _ Mr2)/, ^ /■ /"^ f^ ^(^PN(j3)f(2), ^^ (2) 



ON N 



PNiri)^fp^{c^N\^2,rs)-4\r,,r,))dv,dv,), (17) 



(2) 

where in analogy with the GCE we have introduced the CE direct correlation function c\r^ as 



with Tid given by (|l5|). Eqs. ([l6[)-([l8|) are the main results of the present work. These equations show that it is 
possible to derive the CE counterpart of the direct correlation function and the OZ equation. The new OZ equation 
is consistent with the constraints (|5|) and (|g) as one can readily show by noting that, independently of the functional 
form of c)y , these constraints lead to identities when used in ( [l6| ) and ( |l7| ) (after an appropriate integration) . 

(2) 

Given an explicit form for the intrinsic free-energy functional J-, one can obtain c)^ via functional differentiation 
and then dpN{i^2)/dN is obtained from (p7|). Finally, using (ff6|) one obtains the total correlation function Hn- In 
general, solving these equations is not an easy task and one must resort to self-consistent procedures. There are some 
cases in which one can obtain the required information by means of simple methods. The simplest situation arises for 
the ideal gas. In this case c)^' = 0, Eq. (17) implies dpN{r)/dN = pN(r)/N and, from ( p^ ) we obtain the well-known 



result /lAT — —l/N for the total correlation function of the ideal gas. 

Another simple situation arises in the uniform limit. Considering that in this limit pn{^) ^ pN = N/V {V is the 
volume of the system), one has dpN /dN — pn/N and thus Eq. ( p7| ) becomes 

hN{ri2) = c'^\ri2) + PN c^Vi3)'iw(?-32)dr3 - -^ (19) 

where |rj — r^l = r^ and we have taken into account that the uniform fluid is both translationally and rotationally 
invariant. This equation should be compared with the GCE-OZ equation (|]), that, in the uniform limit, p(r) -^ p = 
Pn, reduces to 

/i(ri2)=c(2)(ri2)+p f c<'^\ru)h{r32)dr3 . (20) 



Then, for uniform fluids, both equations only differ by the 1/A'^ term which vanishes in the thermodynamic limit. 

It is well-known [|ll|,|2| that the GCE radial distribution function g{r) = h{r) + l can be related to gN{r) — /iAr(r) + l 
by means of a series expansion in powers of 1/A^. Using this expansion one has for the total correlation function (with 
PN = p) 

where S{k) = 1 + ph{k) is the structure factor, being h{k) = J h{r) exp(ik • r)(ir the Fourier transform of h{r). In 
what follows we shall indicate how one can rederive this result by means of the OZ equations ( [ig| ) and (|2^) and an 
approximate expression for c)q in terms of c'^'^\r) (conversely, Eq. (pi|) and Eqs. (n9h-(E0) can be used to obtain an 



approximate expression for Cp/)- Since the OZ equations (19) and (|20|) involve convolutions it is convenient to work 



in Fourier space. Subtracting ( pO| ) from ([l9|), after some algebra, we obtain the exact relation 

~ -87r^N-^Sik)S{k) + Acik)SHk) .„„. 

^^(') = 1 - pSik)Ac{k) ('') 

where Ah{k) = h^lk) ~ h{k) and Ac(fc) = cjy (fc) — c(^)(fc). From this equation [or, alternatively, from Eq. (0)] one 
easily obtains Ah{0) — ~S{0)/p, i.e., although Ah{r) is of order l/N, its integral over the total volume V does not 
vanish with increasing N. A simple analysis of Eq. (E2|) can be done by considering that, since the difference between 
the CE and the GCE free-energy functional is of order logiV |jl^, one has that A(P'\k) is of order 1/A^, and thus, 
to first order in AcP'^k) we obtain 

Am « -5!i!«M + (^ _ ^^^sik)m\ ^^(,)^,(,) _ (23) 

We note that, to zeroth order, this equation reduces to Ah{k) w —8Tr^S{k)6{k)/N and thus Ah{r) w ~S{0)/N, i.e., 
the first contribution to Eq. (QW . In order to evaluate (ESh we need an expression (valid to first order) for Ac*^^^ , this 
expression can be obtained from the above-mentioned saddle point approximation for the CE free-energy functional 
fR3 which leads to 



From this approximation, a rather lengthy calculation shows that Eq. ( |23| ) becomes the series expansion result of 
Eq. (pi|). This demonstration is not only a proof of the usefulness of the CE-OZ equation but also it gives further 
support to the saddle point approximation of Rcf. [|0| . We note that this approximation has been shown [Q to yield 
results for the inhomogeneous density that are equivalent to first order to those of the series expansion of the CE 
inhomogeneous density introduced by Gonzalez et al |14[| . 

In summary, we have seen that the fixed- iV constraint of the CE is inconsistent with the standard OZ relation in 
the GCE. By resorting to DFT methods we have identified the source of this inconsistency with the lack of a strict 
one-to-one relation between the density and the external potential which leads to a singular Hessian matrix. This 
problem has been solved by means of an extended variable-space in which we rederive the usual CE Euler-Lagrange 
equation with the advantage that the new Hessian matrix is no longer singular and allows to introduce a CE-OZ 
relation. For both the ideal gas and the uniform fluid this CE-OZ relation only differs from GCE-OZ by a 1/A^ term. 
The presence of this term leads to the exact total correlation function for the ideal gas. For the uniform fluid we 
recover the result of the standard series expansion of the CE total correlation function. 
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